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1.  Introduction 


1.1  The  Nature  of  Ihe  Problem 


Large  regions  of  the  ocean  (amounting  to  about  30%  of  the  earth's  surface) 
are  inadequately  surveyed  for  gravity  data  (Figure  1).  To  improve  present  models 
of  the  global  geoid  it  is  necessary  to  rectify  this  situation  and  to  extend  our  know- 
ledge of  the  gravity  field  into  these  regions,  A number  of  different  methods  have 
been  suggested  for  this  purpose,  including: 

(i)  least-squares  prediction  using  the  auto-covariance  function  of  the 
known  global  gravity  field.  This  is  basically  an  extrapolation 
process,  with  all  the  inherent  shortcomings  of  extrapolation.  Some 
researchers  have  preferred  to  adopt  a zero  anomaly  value  with  a 
large  variance,  rather  than  use  a predicted  value  (Gaposhkin,  1973) 
so  that  the  estimated  value  has  little  influence  on  the  determination 
of  the  potential  coefficients; 

(ii)  computation  of  Ag  from  the  geoid  undulation  (or  its  approximation) 
measured  by  altimetry  from  GEOS-3,  using  either  collocation 
(Rapp,  1977)  or  by  a solution  of  the  inverse  Stokes'  Equation 
combined  with  collocation  (Gopalapillai,  1974); 

(ili)  methods  which  imply  gravity  information  from  known  geophysical 
properties  of  the  ocean  floor.  The  correlation  of  these  properties 
with  gravity  anomalies  can  be  derived  from  areas  of  good  gravity 
data,  and  this  correlation  is  applied  to  known  geophysical  data  in 
areas  lacking  in  gravity  information.  Thus  an  improved  estimation 
of  the  gravity  field  in  this  area  is  obtained. 

It  is  the  purpose  of  this  report  to  concentrate  on  this  last  method.  After 
all,  in  this  approach  one  is  looking  for  a direct  cause  and  effect  relationship,  and 
if  it  is  found  that  the  anomalous  gravity  field  is  closely  related  to  some  widely 
known  geophysical  property,  such  as  the  depth  of  the  ocean,  then  this  information 
can  help  greatly  in  the  extension  of  the  gravity  field  into  unsurveyed  areas. 


1.2  Requirements  for  a Successful  Prediction 


Woollard  and  Strange  (1966,  p.  96)  have  succinctly  outlined  the  elements 
needed  for  a successful  prediction.  These  are  repeated  here. 

"The  degree  of  success  achieved  in  any  prediction  process  depends  on 
the  following; 
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Present  Status  of  Worldwide  Gravity  Coverage 


(a)  liecognition  of  factors  influencing  the  quantity  to  be  predicted. 


(b)  Recognition  of  the  magnitude  of  the  effect  ol  each  individual  factor  1 

on  the  quantity  to  be  predicted.  | 

(c)  Recognition  of  the  interaction  between  individual  factors  and  know-  j 

ledge  of  the  degree  of  interdependence  between  factors. 

(d)  Availability  of  data  concerning  each  factor. 

(e)  One's  philosophy  of  approach  as  to  whether  a given  factor  varies  in 
a random  or  a discrete  way.  " 

It  is  proposed  to  review  the  literature  which  has  explored  the  correlations  ] 

between  various  geophysical  phenomena  in  order  to  gain  some  insight  into  the  i 

various  factors  which  have  been  used  as  correlators  with  the  gnivity  anomaly.  \ 

Much  of  this  literature  is  concerned  with  geophysical  inter}) rotation,  i.e.  the 

inference  of  geophysical  properties  of  the  ocean  floor  and  its  sub-straUi  from  | 

the  gravity  anomalies.  This  inference  is  in  the  op})ositc  direction  to  the  j 

prediction  of  gravity  from  geophysical  })henomena  being  treated  in  this  report.  * 

However,  the  two  approaches  are  very  closely  connected  in  that  techniciues  of  i 

gravity  measurement  and  of  the  modelling  of  the  geophysical  properties  are  common  \ 

to  both,  and  what  is  learned  from  the  geophysical  research  can  certainly  be  used  in  i 

the  gravimetric  applications.  \ 

i 

Before  reviewing  the  literature,  it  is  necessary  to  define  the  terms  and  ; 

explain  the  models  which  are  used  in  the  discussions  of  ocean  floor  geophysics  ' 

that  follow. 


1.3  Definitions 


a)  The  Earth's  Crust  and  Upper  Mantle 

A simple  model  of  the  earth's  crust  and  uppjer  mantle  is  illustrated  in 
Figure  2.  It  will  be  noticed  that  three  principle  zones  can  be  identified;  a sedi- 
mentary layer,  the  crust  consisting  of  both  granitic  and  basaltic  material  and  the 
up)per  mantle  of  ultrabasic  rock.  The  sedimentary  layer  in  mid-ocean  areas  varies  in 
depth  (0. 1 km  to  1 km)  and  density  and  in  some  areas  is  non-existent.  The  mean 
density  of  the  upper  crust  is  usually  taken  to  be  2.  67  gm  cm~^  (cf.  2.  7 to  2.9  gm  cm”‘ 
in  Figure  21.  Its  thickness  decreases  with  increase  in  distance  from  the  continenfcil 
margin,  and  the  ocean  basin  floor  is  found  to  be  between  5 to  10  km.  The  upper 
mantle  is  usually  assumed  to  have  a density  of  3.2  to  3.3  gm  cm"®  (Pick,  Picha 
and  Vyskocil,  1973,  p.  212-214;  Vogt,  Schneider  and  Johnson,  1969,  p.  557). 
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SCHEMATIC  VIEWS  OF  CONTINENTAL  AND  OCEANIC  CRUSTS 


Figure  2:  Schematic  View  of  Continental  and  Oceanic  Type  of  Earth 
Crust  (from  Pick,  Plcha  and  Vyskocil,  1973,  p.  213) 
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It  is  necessary  to  explore  the  theories  of  the  tector’c  processes 
of  the  ocean  floor  in  order  to  appreciate  the  inter-relationships  between  the 
various  geophysical  parameters  which  will  be  referred  to  below  (Section  2.1). 
Since  the  early  1960's  (Hess,  1960;  Vine  and  Mathews,  1963),  it  has  been 
thought  that  new  material  forming  the  crust  may  be  generated  along  the  sub- 
oceanic  ridge  lines  (see  Figure  3).  This  theory  is  biised  partly  on  the  evidence 
provided  by  the  changing  magnetic  orienbition  of  crustal  material.  "The  ocean 
crust  has  been  created  by  dike  injections, . .at  the  axis  of  the  mid-oceanic 
ridges.  As  the  mantle  derivatives  are  injected  into  the  axial  crust  and  cool, 
they  mustacc;uire  a magnetization  depending  upon  the  ambient  strength  and 
direction  of  the  geomagnetic  field..  .The  order  of  magnitude  of  sea-floor  motion 
had  to  be  several  cm/year  if  continental  movements  were  to  be  explained  by 
sea  floor  spreading.  " (See  Vogt,  Schneider  and  Johnson,  1969,  p.  557.) 

An  associated  phenonenon  is  that,  as  the  crust  spreads  from  the  axial 
zone  of  crustal  genesis  and  cools,  it  subsides  (see  Figure  “?).  Hence,  as  will  be 
seen  later  (Section  2.2  (iv))  a relationship  is  found  to  exist  between  crustal  agc' 
and  ocean  depth  thus  establishing  one  of  the  interactions  mentioned  in  (c)  of 
Section  1.1. 


(b)  Accuracy  of  Gravity  at  Sea 

Early  gravity  measurements  at  sea  had  an  expected  accuracy  of  ±5  to 
±8  mgal  (Anderson,  1962,  p.  54),  i.e.  about  one  order  of  magnitude  lower  than 
for  land  measurements.  This  was  mainly  due  to  errors  in  position,  the  motion 
of  the  ship  bearing  the  gravimeter,  the  Eotvils  effect  and  the  difficulty  in  esti- 
mating the  drift  of  the  gravimeter  accurately  (Nettleton,  1976,  p.  116-124). 

More  recently,  the  use  of  the  gyro-stabilized  platform  and  satellite  navigation 
has  improved  this  accuracy  to  abcxit  ±3  mgal  (Lucas,  1971,  p.  H-9;  Khan  et  al. , 
1971,  p.  1). 

Observed  gravity  is  corrected  for  the  usual  latitude  and  drift  factors,  and 
also  a (potentially)large  Eotvds  correction  must  be  applied.  The  result  is  the  free- 
air  anomaly  as  gravity  is  measured  at  sea  level  (effectively).  Individual  measure- 
ments can  then  be  combined  in  some  standard  technique  to  obtain  mean  free-air 
anomalies  for,  say,  the  5°x  5°  or  l°x  1°  blocks  which  are  required  for  global 
gravity  data  banks.  The  accuracy  of  these  mean  values  is  dependent  on  the  nature 
of  the  ocean  floor,  the  amount  of  data  and  the  method  used  in  obtaining  the  mean, 
but  for  a well  represented  l°x  1°  block  the  accuracy  can  be  less  than  i 10  mgals. 
For  areas  with  little  or  no  known  data,  an  accuracy  of  about  three  times  this  is 
possible;  Rapp  and  Rummel  (1975)  assumed  a sUindard  deviation  of  ± 30  mgals 
for  l°x  1°  means  in  the  U.  S.  Calibration  Area,  this  Ix^ing  the  r.  m.s.  of  the 
free-air  gravity  anomalies  in  the  area. 


(c)  Bouguer  and  Isostatic  Anomalies 
(i)  Bouguei-  Anomalies 


The-  normal  Bouguer  correction  implies  the  removal  of  an  infinite  slab 
of  material  of  assumed  density  whose  thickness  is  the  height  of  the  observation 
point  above  mean  sea  level.  This  correction  enables  one  to  compare  gravity 
anomalies  at  sea  level  after  the  attraction  of  material  above  this  level  has  been 
removed.  The  aim  of  a 'Bouguer'  correction  at  sea  is  similar.  The  intent  is 
to  remove  the  effect  of  the  density  contrast  between  seawater  and  bedrock  so 
that  gravity  effects  from  sources  below  the  sea  floor  are  more  apparent. 

The  correction  is  computed  in  the  usual  way  except  that  the  density 
contrast  between  sea -water  and  rock  (usually  taken  to  be  1.  64  gm  cm  is 
used  (Nettleton,  1976,  p.  280).  Thus; 


(1.1)  Ag0)uo  - AgFA  + 0.0687.  h' 

where  h'  is  the  depth  of  the  ocean  (see  Figure  3).  To  compute  the  l°x  1°  mean 
value  the  mean  free-air  anomaly  and  mean  depth  are  used  instead  of  the  discrete 
value.  Terrain  corrections  could  also  be  applied  to  correct  for  the  fact  that 
adjoining  l"  x 1°  depths  are  different  to  the  computation  value  and  thus  the  model 
of  an  infinite'  slab  used  in  computing  the  correction  is  inadequate. 


(ii)  Isostatic  Anomalies 


The  theory  of  isostacy  assumes  that  there  is  a state  of  equilibrium  at  a 
certain  level  (T  = 30  km)  below  the  earth's  surface.  In  the  Airy  version  of  the 
theory  this  is  achieved  by  variations  of  thickness  in  the  crust,  which  is  assumed 
to  have  a uniform  density  throughout.  Thus  in  areas  of  low  relief  (A  in  Figure  3) 
the  crust  is  thinner  than  normal,  and  vice  versa  in  areas  of  high  relief  (B  in 
Figure  3),  The  net  effect  is  that  hydrostatic  pressure  at  depth  T is  equal  under 
both  A and  B (see  Heiskanen  and  Vening  Meinesz,  1958,  p.  135-142). 

The  thickness  of  the  anti-root  (f  in  Figure  3)  iS; 

■-I 

« 

(1.2)  t'  = (Pc  - /0.t)/Aph' 

where 

• ^ 

Pc  ^ density  of  the  crust  (^  2. 67  g cm~^) 

Pm  = density  of  sea  water  (^  1.03  gcm“^) 

i 
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Ap  = density  contrast  between  crust  and  mantle  (=-0.6  g cnT^) 
h'  = depth  of  water 

The  total  thickness  of  the  crust  under  the  ocean  (To)  is  therefore: 


(1.3) 


To  = T - t'  - h' 


Corrections  can  now  be  applied  to  the  observed  free-air  gravity  anomaly 
to  account  for  the  fact  that  the  column  beneath  the  observation  point  differs  from 
a standard  column  of  30  km  depth  and  p=  2.  67  g cm”^. 

Other  theories  using  the  theory  of  isostasy  will  not  be  treated  here.  The 
interested  reader  is  referred  to  (Ibid.,  p.  131-142)  and  Kivioja  (1963,  p.  17-22). 


2.  Geophysical  Factors  Related  to  Gravity  Anomalies  in  Ocean  Areas 
2.1  Introduction 


The  oceans  are  composed  of  three  morphological  units:  continental 
margins,  mid-oceanic  ridges  and  ocean  basins.  Mid-oceanic  ridges  constitute 
about  30%  of  the  oceanic  area  (Vogt,  Schneider  and  Johnson,  1969,  p.  362)  and, 
because  of  their  role  in  the  tectonic  process  (Section  1.3(a)),  present  a complex 
picture  in  the  interpretation  of  anomaly -causing  structures.  The  ocean  basins 
may  possess  features  such  as  abyssal  plains,  mid-ocean  canyons,  ocean  rises, 
micro -continents,  seamounts  and  fracture  zones  (Ibid.,  p.  562-573).  All  but 
the  last  feature  are  either  tectonically  quiet  or,  as  in  the  case  of  seamounts, 
are  isolated  and  may  have  little  effect  on  mean  values  of  height  and  gravity  for 
l°x  1°  blocks.  The  fracture  zones  can  have  a depth  of  1 km  and  may  cover 
extensive  areas.  They  involve  crust  and  upper  mantle,  and  depending  on 
magnitude  and  extent,  are  likely  to  complicate  the  relationship  which  may 
otherwise  exist  between  geophysical  features  of  the  ocean  floor  and  the  gravity 
field. 


It  is  necessary  to  place  limits  on  the  extent  of  the  present  investigation 
into  the  correlation  between  gravity  anomalies  and  geophysical  properties.  Kor 
example,  there  seems  little  value  in  reviewing  findings  at  continental  margins  as 
most  of  these  are  already  well  surveyed  for  gravity  (as  witnessed  by  the  fact  that 
the  literature  deals  heavily  with  geophysical  interpretation  from  gravity  in  these 
regions),  investigation  into  the  fine  structure  in  localized  areas  (such  as  the 
Hawaiian  Emperor  Seamount  chain,  see  Lucas,  1971;  Watts,  1976)  will  imobably 
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have  limited  application  and  will  generally  not  be  included  in  the  survey.  Papers 
which  explore  the  gravity  relationships  to  such  tectonic  properties  as  lithospheric 
flexure,  plate  thickness,  mantle  convection  and  heat  flow  are  also  omitted  as 
data  on  these  matters  is  limited  to  areas  well  surveyed  for  gravity,  and  the 
conclusions  concerning  their  relationships  are  still  a matter  of  conjecture. 

The  obvious  place  to  concentrate  research,  then,  is  in  ocean  basins 
which  arc  uncomplicated  by  fracture  zones  or  large  sea  mounts.  In  the  section 
which  follows,  a number  of  papers  are  reviewed.  These  are  not  limited  to  ocean 
basins  but  are  broader  in  scope  so  that  general  relationships  (and  their  compli- 
cations) can  be  seen  in  perspective.  Following  this  review  some  numerical  tests 
will  be  made  in  an  ocean  basin  area  to  help  evaluate  some  of  the  relationships 
which  are  found  to  exist  in  the  literature. 


2.2  Review  of  Relationships  Between  Gravity  Anomalies  and  Geophysical 
Phenomena 


(a)  Global  Relationships 

A number  of  researchers  have  analyzed  the  correlation  between 
geophysical  or  geomorphological  phenomona  and  the  global  gravity  field  in  ocean 
areas  of  the  earth.  This  is  a broader  approach  than  is  required  in  this  present 
investigation  but  may  be  useful  in  providing  a 'bias'  (or  datum  shift)  to  predicted 
anomalies  in  ocean  areas  devoid  of  gravity  data. 

Kivioja  and  Ijewis  (1966)  used  worldwide  5°x  5°  mean  elevations  and 
depths  in  an  Airy-Heiskanen  isostatic  model  to  compute  free-air  anomalies  for 
each  5°x  5”  block.  Although  no  comparison  with  known  gravity  data  is  given, 
the  geoid  majis  computed  by  Stokes  formula  from  both  the  generated  and  known 
data  were  compared.  The  authors  concluded  that  a large  part  of  geoid  undulation 
is  due  to  bathymetric,  topographic  and  isostatic  masses,  and  felt  that  the 
differences  n-vealed  information  about  the  hidden  density  distributions  inside  the 
earth.  It  is  dangerous  to  draw  too  broad  a conclusion  from  this  study,  however, 
as  both  solutions  contained  common  data,  the  generated  data  being  used  in  areas 
where  there  was  no  surveyed  data. 

* 

Wool  lard  (1969,  p.  286)  noted  that  ge^idal  anomalies  defined  by  the 
then  current  harmonic  representation  of  satellite  data  bore  no  relationship  to 
surface  mass  distribution,  but  notes  a correlation  between  surface  gra.  ity  and 
tectonic  activity.  * 

Ijambeck  (1971)  also  pointed  out  the  presence  of  positive  free-air  anomalies 
over  19  ridges  selected  from  all  major  oceans  of  the  world,  and  developed 
lithospheric  models  to  explain  the  variation  in  gravity  over  these  ridges  (Ibid., 


-8- 


p.  43-48).  Anderson  et  al.  (1973)  confirmed  this  correlation,  but  found  no 
correlation  between  spreading  rates  and  observed  crusUil  depths.  Roufosse 
and  CJaposhkin  (1976)  also  noted  the  correlation  of  Ag  with  mid-oceanic  ridges  and 
also  in  volcanic  areas  but,  when  comparing  short  wavelength  features  of  the  global 
Ag  field  with  topography  concludes  that  the  comdation  was  jKior,  particularly  in 
ocean  regions  (see  also  Roufosse,  1977). 

It  is  obvious  from  this  survey  that  there  is  no  simple  general  relationship 
which  could  be  used  on  a global  scale  to  predict  Ag  from  geophysical  data  in  ocean 
areas.  As  recognized  by  Khan  et  al.  (1971),  it  is  necessary  to  restrict  the 
study  to  tectonically  quiet  and  uniform  areas  in  order  to  limit  the  varial)les  which 
influence  the  gravity  values  in  (x:ean  areas. 


(b)  Relationship  with  Topography  or  Bathymetry 

Many  investigators  have  noted  a correlation  between  depth  of  water  (or 
sea  floor  topography)  and  gravity  in  various  ocean  regions.  These  relationships 
are  complex  when  over  tectonically  active  ridges.  But  in  ocean  basins  one  would 
expect  positive  gravity  anomalies  in  shallower  than  normal  areas  (and  vice  versa) 
providing  the  density  of  the  sub-strata  rcmiains  more  or  less  uniform  throughout 
the  region. 

Taiwan!  and  Le  Pichon  (1969)  analyzed  5° x .'i”  mean  anomalies  and  depths 
in  both  the  North  and  South  Atlantic  and  found  a very  strong  correlation  l)ctween 
topography  and  gravity  in  the  North  Atlantic  while  in  the  South  Atlantic,  the 
correspondence  was  less  apparent.  (The  gravity  field  in  the  South  Atlantic  at  that 
time  was  not  well  defined.)  Mathews  et  al.  (1969)  in  a more  localized  analysis  in 
the  Peak  and  Freen  Deeps  in  the  North  Atlantic  felt  that  "large  anomalies. . .are 
mainlj  due  to  the  topography  and  that  there  is  an  almost  uniform  mass  distribution 
below  the  sea  floor"  (Ibid.,  p.  533).  Dehlinger  (1970)  corroborates  this  relation- 
ship in  the  North  Pacific  ocean  basins,  noting  that  most  of  the  topographic  highs 
(usually  seamounts)  are  characterized  by  fx)sitive  free-air  anomalies,  although 
he  felt  that  some  of  the  Ag's  eould  not  be  correlated  with  topographic  features, 
and  that  the  correlation  was  more  on  a local  than  a regional  l)asis.  In  fact,  he 
concludes  that  topography  is  a poor  to  unsatisfactory  guide  for  estimating  regional 
anomalies,  even  where  the  extent  of  isostatic  eejuilibrium  has  been  determined 
(Ibid.,  p.  36.3). 

It  has  been  suggested  that  this  relationship  between  ocean  floor  topography 
and  gravity  anomalies  holds  even  in  tectonically  active  regions.  F^arly  studies  over 
the  mid-oceanic  rise  of  the  North  Atlantic  showed  that  isostatic  equilibrium  existed 
in  this  ridge  area  and  that  topography  accounted  for  most  of  the  anomalies  here 
(van  Andel  and  Bowen,  1968).  And  in  1973,  Anderson  and  others  showed  a positive 
correlation  between  free-air  anomalies  and  differences  in  crusUil  depths  of  the 
mid-oceanic  ridge  systems.  However,  they  found  no  correlation  between  spreading 
rate  and  gravity,  and  felt  that  no  uniform  txdationship  held  for  all  (x;eans  between 
spreading  rate  and  observed  crustal  depths  (Anderson  et  al. , 197.3). 
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The  fort'KoinK  suKSests  there  is  correlation  between  sea-floor 
to|>onraphv  and  j^ravity  anomalies,  jxirticularly  in  the  Atlantic  Ocean  and  for 
heal  areas  s^nerally.  However,  it  appears  that  no  simple  relationship  holds 
for  la  rue  areas  in  the  I’acifie  Oc«‘an. 


(c)  Kelationship  with  Depth  Anomalies 

Th<'  depth  anomaly  is  defined  as  the  difference  between  the  measured 
depth  and  a depth  value  neneratt'd  by  an  empirical  model  based  on  crustal  age 
(Anderson  <•(  al. , 197.1,  p.  4().‘{).  Sclater  et  al.  (1975)  examined  the  correlation 
between  l°x  1°  and  .5'^x  .5°  mean  depth  anomalies  and  gravity  anomalie.s  in  the 
relationship  iK'tween  these  two  data  sets.  However,  for  the  l°x  l"  means,  the 
correlation  aveniged  thnxigh  the  whole  area  was  poor  and  could  not  be  used  with 
confidence  in  prediction  (Ibid. , p.  10.16),  although  certa in  local  a reas  (e.  g.  South 
of  Iceland  and  near  the  Azores)  show  strong  correlation.  Marsh  and  Marsh  (1976) 
noticed  visual  correlations  between  the  residual  depth  anomalies  and  Ag  in  the 
central  and  eastern  Pacific  oceans.  They  also  found  a series  of  linear  positive 
and  negative  'residual'  anomalies  across  the  Pacific.  (A  residual  anomaly  was 
defined  here  as  that  field  remaining  after  a 12  degree  and  order  field  is  sub- 
tracted from  a field  defined  to  degree  and  order  22).  This  suggested  to  them 
the  pre.sence  of  longitudinal  convective  n)lls  beneath  the  Pacific  plates  which  may 
be  significant  when  developing  a prediction  mcxiel  for  the  Pacific  region. 


(d)  Helationships  with  Crustal  Age 

The  theory  of  sea  floor  spreading  provides  a model  which  enables  one  to 
understand  the  interrelationships  between  some  geophysical  parameters.  As  shown 
in  Sc‘cli«)n  1. 3 (a)  both  the  magnetic  profierties  and  depth  of  the  ocean  floor  can  be 
seen  as  a tanction  of  age.  Hoth  these  properties  are  relatively  easy  to  measure 
and  a large  body  of  data  is  available  for  them.  l‘\irthermore,  as  reported  in 
Woollard  and  Daugherty  (1974,  p.  5),  Woollard  has  found  "systematic  interrela- 
tionships between  depth  of  water,  thickness  of  crystalline  rfjck  crust,  mean 
velocity  of  the  crust  and  velocity  of  the  mantle  in  the  Pacific  region,  and  has 
observed  that  these  interrelationships  change  with  crustal  age  and  past  crustal 
spreading  rate". 

In  one  of  the  few  attempts  to  predict  gravity  from  other  geophysical 
information,  Daugherty  has  used  crustal  age  derived  from  (i)magnetism  and 
(ii)  water  depth  to  predict  l”x  1°  and  5°  x 5°  mean  gravity  anomalies  in  the 
northeast  Pacific.  The  1°  x 1"  areas  involved  are  classified  according  to 
U)|K)graphic  or  tectonic  type  in  (Daugherty,  1975,  p.  8),  and  shows  that  about 
6.5'iri  of  these  areas  are  undisturbed  by  major  topographic  or  tectonic  features. 

His  prediction  method  is  outlined  below. 
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(i)  An  area  in  the  north  Pacific  bounded  roujjhiy  by  6°  - 1°, 

210°  ■ X ^ 255°  was  chosen  as  a test  area  as  it  was  c<»nsidered  to  have  adei)uaU‘ 
data  and  had  relatively  uncomplicated  topographic  and  tectonic  characteristics, 

(ii)  The  1°  x 1°  mean  free~air  gravity  anomalies  (^f*)  wen  obtained 
from  the  Department  of  Defence  (DoD)  gravity  library.  (Xiean  floor  aj  cs  an<l 

were  mapped,  and  the  Agf*  's  averaged  over  5 million  year  (m.\  .)  inU-rvals. 

(iii)  These  5 m.y.  averages  were  now  plotted  against  age  up  to  80  m.y., 
beyond  which  lack  of  data  precluded  any  meaningful  analysis.  A serie:.  of  poly- 
nomials were  now  fitted  to  this  plot  and  the  best  fitting  polynomial  assumed  to 
represent  the  Ag  vs.  age  relationship. 

(iv)  l°x  1°  mean  values  of  AgfA  wert?  now  'predicted'  from  the  l^’x  1° 
mean  crustal  ages.  A comparison  of  the  predicted  vs.  known  Ag 's  en  ibled  a 
statistical  assessment  of  the  success  of  the  prediction. 


Itesults  and  Comments 


The  root  mean  square  (RMS)  of  the  differences  in  (iv)  ranged  rom  ±10.5 
mgal  to  ± 12  mgal,  the  former  for  when  ages  were  interpolated  directh  from  an 
age  map  and  the  latter  computed  from  the  age-depth  relationship  referi  c*d  to  in 
(c)  above.  This  is  abixit  a 40%  improvement  over  the  RMS  of  the  graviiy  anomalies 
themselves  (±17  mgal)  and  compares  favorably  with  the  accuracy  of  th«  'known' 
mean  anomalies,  which  varies  between  ±5  to  ±25  mgal  and  has  a mean  of  about 
±13  mgal. 

Unfortunately,  the  success  of  this  technique  was  limited.  When  applied 
to  a new  data  set  in  the  South  East  Pacific  (where  nearly  50%  of  the  set  had 
uncertain  topographic  type)  it  gave  disappointing  results,  and  it  was  ne'-essary  to 
conclude  that  in  this  region,  there  is  no  apparent  correlation  between  i\g  and  age 
(Daugherty,  1975,  p.  19).  The  complexity  of  the  ocean  floor  structure  (Woollard 
et  al. , 1975,  p.  3)  and  the  uncertainty  of  the  gravity  data  (Daugherty,  1975,  p.  19) 
were  cited  as  possible  reasons  for  this.  In  fact,  the  points  (a)  to  (d)  in  Section  1.2 
remained  unsatisfied  and  an  unsuccessful  prediction  resulted. 


2.3  Conclusions 

It  appears  that  there  is  a good  correlation  between  depth  and  free-air 
gravity  anomalies  in  ocean  basins.  A more  fundamental  relationship  bt  tween 
anomalies  and  crustal  age  has  been  suggested  and  has  been  used  successfully  in 
one  ocean  area.  However,  there  appears  to  be  no  direct  relationship  which  holds 
generally  in  all  ocean  areas,  and  the  inadequacy  and  inaccuracy  of  the  < bservixl 
gravity  itself  precludes  further  investigation  in  this  field.  The  unpromising  nature 
of  these  comments  should  be  viewed  in  the  light  of  Section  4 where  reference  is  made 
to  alternative  methods  of  predicting  anomalies  in  ocean  areas. 
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C«)lh>cation  of  Gravity  Anomalies  from  Crustal  Age 


1 Introduction 


Prediction  techniques  using  collocation  in  either  the  univariate  or  multi- 
variaU'  mode  (sec  Moritz,  1972)  are  now  widely  used  in  predicting  potential-related 
parameters.  However,  in  practice,  their  use  has  been  limited  to  parameters  such 
as  anomalies,  geoidal  undulations  or  deflections  of  the  vertical,  which  are  directly 
related  to  anomalous  earth  potential.  The  question  now  is  whether  or  not  it  is 
possible  and  feasible  to  predict  gravity  by  collocation  from  some  geophysical 
parameter  which  is  not  directly  related  to  the  anomalous  potential. 

The  theoretical  justification  for  this  can  be  found  in  Tscheming  (1974), 
where  relations  assigning  a density  distribution  to  the  harmonic  part  of  the  earth 
are  derived.  Tscheming  proves  it  is  possible  to  combine  the  values  of  density 
anomalies,  gravity  anomalies  or  other  gravimetric  quantities  in  least  squares 
collocation  using  the  covariance  models  derived  therefrom  (Ibid.,  p.  13-14). 

This  derivation  is  based  on  a global  model,  but  it  is  not  unreasonable  to 
apply  the  technique  to  local  fields  also.  The  results  of  such  an  application  are 
described  below. 


3.2  Prediction  of  Free-Air  Gravity  Anomalies  by  Collocation 
(a)  Covariance  Analysis 

l-’or  purposes  of  comparison  the  data  set  used  by  Daugherty  (1974, 
Appendices  I and  II)  were  used  for  analysis.  The  auto-covariance  functions 
(a  c f)  of  the  free-air  anomaly  (^rA)  and  of  crustal  age  (CA)  and  the  cross- 
covariance function  (cef)  of  Ag^A  with  CA  were  computed  from  the  482  l°x  1° 
areas  means  and  the  resulting  functions  illustrated  in  Figures  4 to  8.  The  units 
of  CA  are  million  years  (m.y.).  It  will  be  noted  that  the  acf  for  behaves 
typically,  from  (()=  0®  to  it  = 1.5®,  but  thereafter  is  very  flat,  reflecting  the 
smoothness  of  the  field.  C (it)  = 0 at  ((i*  15? 5 which  again  shows  that  there 
exists  positi\e  correlation  between  anomalies  over  much  larger  separations 
than  normall  >'  exist  over  land  areas.  The  acf  for  CA  falls  almost  linearly  from 
a C (0)  value  of  400  m.y.’  and  attains  zero  at  (()=>•  16.  5°.  The  cef  of  AgrA  with 
CA  is  also  flat,  crossing  C {4>)  at  i/)^  15? 5. 

The  acf  for  ocean  depth  is  shown  in  Figure  7 and  is  similar  to  that  for 
CA,  at  least  to  10°,  after  which  it  increases  slightly  and  then  continues  to 
decompose  slowly.  The  cef  of  Agr,  with  depth  is  atypical,  reaching  its  max- 
imum at  S°  after  a local  minimum  at  2°.  Such  behaviour  is  not 
uncommon  in  cef's  and  was  noted  by  this  author  in  the  cef  of  AKfa  with  geoidal 
undulations  in  the  U.  S.  Calibration  Area  (Kearsley,  1977).  It  is  surprising  to 
find  the  lag  at  which  this  maximum  occurs  is  so  large.  It  is  hard  to  accept  the 
implication  f i-om  this  that  values  8®  from  a computation  point  have  effectively  the 
same  influence  on  the  prediction  as  values  adjacent  to  it.  Use  of  such  functions 
should  be  tempered  with  critical  caution. 
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(b)  Prediction  Results 


The  prediction  technique  is  given  in  Moritz  (1972)  and  will  not  be 
described  here.  In  the  prediction  the  5 closest  l®x  1,°  mean  crustal  ages  (as 
interpolated  from  the  age  map)  were  used  to  find  the  mean  gravity  and  its  estimated 
error  of  the  subject  l^x  area.  The  difference  between  the  predicted  and  'known' 
AgFA  value  were  used  to  find  accuracy  estimates.  It  should  be  remembered  that 
the  '5  closest'  CA's  would  include  the  CA  for  the  subject  area  itself,  but  the 
remaining  four  may  not  be  adjacent  to  this  area  as  the  data  set  was  incomplete. 

The  results  are  illustrated  in  Figure  9.  The  RMS  of  the  differences  is 
±11.5  mgal,  which  is  not  quite  as  good  as  that  obtained  by  Daugherty  (±10.4  mgal) 
but  is  certainly  of  the  same  order.  It  compares  well  with  the  theoretical  error  of 
prediction  (=>-11.2  mgal).  The  RMS  of  the  anomalies  themselves  is  ±17.2  mgal, 
so  the  prediction  represents  an  improvement  of  about  30%  over  an  assumption  of 
a zero  anomaly  for  the  whole  data  set. 

It  must  also  be  remembered  that  the  accuracy  of  the  'known'  gravity  data 
is  at  times  poor  (up  to  ±25  mgal),  averaging  only  about  ±13  mgal.  An  RMS  of 
differences  of  ±10  to  ±11  mgal  is  therefore  quite  good. 

A further  comparison  was  made  between  'known'  values  and  values 
predicted  from  the  AgpA  values  of  the  5 nearest  l°x  1°  areas,  excluding  the 
subject  area.  The  distribution  of  the  differences  is  shown  in  Figure  10.  The  RMS 
in  this  case  was  ±9.  3 mgal,  i.e.  of  the  same  order  as  the  prediction  from  crustal 
ages.  It  docs  appear  that  prediction  of  Ag  from  crustal  age  is  successful  in  this 
area,  approaching  the  limit  set  by  the  accuracy  of  the  known  gravity  data. 


3.3  Testing  the  Use  of  Bouguer  and  Isostatic  Anomalies 
(i)  Bouguer  Anomalies 

It  is  common  practice  to  use  Bouguer  anomalies  when  predicting  in  local 
areas  over  land,  as  this  field  is  smoother  and  and  more  amenable  to  linear  inter- 
polation. The  computation  of  'Bouguer'  anomalies  at  sea  has  already  been  outlined 
in  Section  1.  •!  and  are  known  to  attain  high  absolute  values  in  ocean  areas  (e.g. 
Heiskanen  and  Vening-Meinesz,  1958,  p.  144,  197).  However,  their  magnitude  is 
unimportant  m this  context.  What  is  needed  for  successful  prediction  is  a smooth 
field  with  sm  ill  residuals  about  the  mean. 

A covariance  analysis  of  the  Bouguer  anomalies  showed  that  these 
anomalies  w<  uld  certainly  not  produce  such  a field.  It  is  obvious  that  the  correction 
values  vary  preatly  from  one  l°x  1°  area  to  the  next  and  in  fact  disturb  an  already 
smooth  field.  The  variance  of  this  field  is  very  large  (2420  mgal®)  and  the  auto- 
covariance function  (acf ) drops  steeply  (cf.  the  free-air  anomaly  acf ) to  cross 
C ((<>)  = 0 at  (t  » 19°.  The  cross-covariance  of  Bouguer  anomaly  with  age  is  also 
large  (about  10  mgal  m.y. ).  It  is  obvious  that  use  of  Bouguer  anomalies  will  pro- 
duce poorer  results  than  those  obtained  using  the  free-air  anomalies,  and  no 
predictions  v ere  attempted  with  them. 
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(ii)  Isostatic  Anomalies 


It  appeared  th;it  isostatic  anomalies  (Akiso)  would  produce  more 
promising  results  than  Bouguer,  and  even  free-air  anomalies,  as  the  isostatic 
correction  results  from  a more  sophisticated  model  which  could  account  for 
some  of  the  'noisy'  features  in  the  data.  The  isostatic  anomalies  for  the  data 
set  were  computed  as  described  in  1 . 3 (c)  and  the  results  subject  to  covariance 
analysis. 


The  isostatic  anomalies  proved  to  be  almost  identical  with  the  free-air 
anomalies,  the  largest  difference  being  =■  1 mgal.  As  a re.sult  the  statistical 
analysis  is  also  practically  identical  (C  (0°)  = 158.3  mgal®  for  Agiao,  cf. 

159.0  mgal*  for  AgfA).  It  was  obvious  that  there  would  be  no  impnivement  in 
accuracy  if  one  used  isostatic  rather  than  free-air  anomalies,  and  no  such 
prediction  was  attempted. 


Comments 

It  appears  that,  in  the  test  area  at  least,  the  ocean  acts  as  a filter' 
to  local  anomalies  in  the  crust,  and  that  the  free-air  anomalies  are  already 
smooth  and  suitable  for  use  in  the  prediction. 

It  also  appears  that  isostatic  equilibrium  holds  in  this  area  as  the 
free-air  and  isostatic  anomalies  are  so  similar. 

Preliminary  calculations  of  the  terrain  correction  showed  that  this  was 
of  such  small  magnitude  that  it  would  not  have  any  impact  in  improving  the 
prediction.  This  is  due  mainly  to  the  size  of  the.  sejja ration  between  adjoining 
blocks  (at  least  110  km)  and  the  depths  of  the  ocean  being  dealt  with  (=-  3 to  5 km). 


4.  Conclusi(ns 

A r<  view  of  the  literature  showed  that  a general  relationship  existed 
between  gravity  anomalies  and  depth  of  water  in  ocean  areas.  Woollard  has 
suggesk'd  that  crustal  age  is  a more  fundamental  parameter  to  use  in  prediction, 
based  on  the  heory  of  sea  floor  spreading.  Early  investigations  using  a poly- 
nomial relationship  between  crustal  age  and  gravity  anomalies  showed  promise, 
giving  an  accuracy  of  ± 10.  5 mgal  for  a l°x  1°  mean  value  (Woollard  and  Daugherty, 
1974).  Howe  er,  extension  of  this  approach  into  a tectonically  uncertain  area 
failed,  highl i!’,hting  the  fact  that  research  in  ocean  areas  is  severely  hampered  by 
poor  coverag.  and  accuracy  of  the  gravity  and  lack  of  knowledge  of  the  nature  of 
the  ocean  flo«  r and  its  sub-strata. 


Collocation  was  used  to  predict  free-air  gravity  anomalies  fnun 
Daugherty's  crustal  age  data.  This  produced  results  comparable  to  th(  test 
mentioned  above  (±11.2  mgal)  and  represented  an  improvement  of  about  30% 
over  an  assumption  of  zero  anomaly  for  this  data  set.  It  also  comfjare  ; 
favorably  with  the  average  accuracy  of  the  known  data  in  this  region  (±  13  mgal). 

(The  uncertainty  of  this  data  reaches  ±23  mgal  in  some  cases,  and  one  must 
question  its  suitability  for  comparison  purposes  in  such  cases.) 

The  collocation  approach  has  several  advantages,  it  is  a sim|)le  matter 
to  extend  the  technique  to  include  a number  of  parameters  (e.g.  gravity  anomalies 
and  crustal  age)  in  the  prediction.  Error  estimates  of  the  data  can  be  incoq)or- 
ated  in  the  solution,  which  also  provides  an  estimate  of  the  accuracy  of  the 
prediction.  Some  preliminary  analysis  of  the  data  is,  of  course,  required  to  find 
the  relevant  covariance  functions. 

Recent  research  has  shown  that  the  need  for  the  type  of  predici  ion  discussed 
in  this  report  has  been  lessened  in  areas  where  reliable  satellite  altimetry  data  will 
be  available.  Developments  in  the  prediction  of  gravity  anomalies  from  f:EOS-3 
altimetry  have  shown  that  this  technique  is  cajiable  of  a superior  accuracy  of  ±8  mgal 
for  l°x  1°  area  means,  and  ±2.5  mgal  for  5°x  5°  means  (Rap)),  1977;  see  also 
Appendix  I).  Many  ocean  areas  of  geodetic  interest  have  been  covered  by  salellite 
altimetry  and,  as  more  data  becomes  available,  will  be  j)rocessed  to  providi'  a 
fairly  complete  ocean  coverage  (except  at  high  latitudes)  to  an  accuracy  approaching 
that  obtainable  from  ship  borne  gravimetry.  The  anomaly  is  deduced  from  direct 
measurement  of  the  geophysical  approach  described  herein.  In  fact,  it  is  probable 
that  in  the  near  future  geoidal  undulations  and  gravity  anomalies  obtained  from 
satellite  altimetry  will  aid  in  the  interpretation  of  large  scale  (-^5°)  or  small  scale 
( -0?3)  geophysical  features  of  the  ocean  floor  and  its  sub-strata  (e.g.  Howin,  1975). 
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Appendix  I - Anomaly  Comparisons 


The  table  below  gives  a comparison  of  a small  sample  of  free -air 
gravity  obtained  from  various  sources  for  l°x  1°  areas  stated. 

Ag'  are  from  the  DoD  library  (DMAAC,  1972,  1972)  used  b\  Daugherty 
and  this  author  for  their  analysis  and  prediction. 

Ag"  are  from  an  updated  version  of  the  DoD  gravity  library  (DMAAC, 
1976). 

Ago  are  values  predicted  by  Daugherty  (1974,  Appendix  II)  from 
crustal  age. 

IV 

AgK  predicted  by  collocation  from  crustal  age  in  this  report  i section  3). 

y 

Ag"  alt  predicted  by  R.  H.  Rapp  (private  communication,  December  1977) 
from  GEOS-3  altimetry. 


X 

'Known'  Data 

A 

Ago 

Predicted  Data 

(Degrees) 

Ag 

A 

Ag 

AgK 

<I 

i 

! 39 

216 

-10 

-20  ±11 

-13 

-11  ±11.2 

- 8 ±9  i 

39 

217 

-17 

-19  ±22 

-16 

-13  ±11.2 

-14  t6  j 

39 

218 

-20 

-12  ±17 

-17 

-13  ±11.2 

-20  ±7  j 

40 

217 

-33 

-51  ±20 

-18 

-14  ±11.2 

-28  ±9  ! 

40 

218 

-16 

-27  ±19 

-19 

-15  ±11.2 

-31  ±8 

40 

226 

-12 

-39  ±13 

-19 

-18  ±11.2 

-32  ±7 

40 

L _J 

233 



-11 

-22  ± 7 

-17 

-20  ±11.2 

-22  ±9 

(all  gravity  values  are  in  mgal). 


It  is  worth  noting  the  large  differences  between  the  'known'  data  and 
the  generally  large  uncertainty  of  the  Ag"  values.  For  (p  = 40°  the  Ag'^AU 
compares  well  with  Ag"  (except  at  X=  217),  while  for  the  remainder  it  compares 
more  favorably  with  Ag' . A comparison  for  such  a small  sample  is  meaningless, 
but  it  does  help  to  point  up  the  uncertainty  of  the  known  data  and  the  apparent 
improvement  in  accuracy  in  the  altimetry -derived  values. 
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